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ABSTRACT 


This  thesis  presents  a  method  for  finding  the  trajectory  to 
complete  an  air-to-air  missile  intercept  which  maximizes  information. 
This  is  accomplished  by  formulating  a  parameter  optimization  problem 
and  using  a  penalty  function -Lagrange  multiplier  method  to  solve  for 
the  optimal  path.  The  performance  index  is  the  trace  of  the  informa¬ 
tion  matrix.  This  information  matrix  is  derived  using  an  extended 
Kalman  filter  formulation,  in  cartesian  .'oordinates,  which  makes  state 
estimates  based  only  on  angle  measurements.  The  trace  operation  on 
the  information  matrix  is  used  because  the  trace  and  the  integration 
operations  commute  allowing  a  scalar  performance  index.  Further,  re¬ 
duction  in  the  functional  form  of  the  performance  index  is  achieved  by 
weighing  the  Information  matrix  by  the  inverse  of  the  measurement  power 
spectral  density.  This  also  avoids  numerical  difficulties  near  inter- 
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SECTION  I 


INTRODUCTION 

Many  of  the  modern,  tactical,  air-to-air  missiles  employ  passive 
seekers  such  as  infrared  heat  seekers  in  their  guidance  systems.  For 
particular  geometries  where  the  line-of-sight  inertial  angle  is  small, 
the  proportional  navigation  system,  which  drives  the  inertial  line-of- 
sight  angle  rate  to  null,  works  reasonably  well.  If  the  intercept 
geometry  is  more  severe,  or  if  the  target  is  accelerating,  then  obtain¬ 
ing  range  and  range-rate  information  becomes  important,  especially  to 
derive  the  time-to-go  to  intercept.  Since  passive  or  jammed  sensors 
may  only  have  angle  information,  range  and  range  rate  can  only  be  ob¬ 
tained  by  extraction  from  sophisticated  filter  design. 

In  Reference  1,  Sammons  et  al  compare  various  implementations 
of  the  extended  Kalman  filter  for  missile  onboard  operations.  This 
filter  uses  angle  measurements  to  estimate  relative  position,  relative 
velocity,  and  target  acceleration.  Estimates  of  range  and  range  rate 
are  physically  obtained  by  having  the  inertial  line-of-sight  angle 
constantly  changing.  However,  the  proportional  navigation  guidance 
law  is  always  attempting  to  drive  the  angle  rate  to  zero.  Since  these 
angles  are  the  only  source  of  information,  it  seems  obvious  that,  when 
the  angle  rate  approaches  zero,  there  will  be  little  or  no  information 
to  the  filter  about  range  and  range  rate.  The  guidance  law  is  working 
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against  the  estimator  by  commanding  a  flight  path  which  denys  informa¬ 
tion  to  the  estimator. 

The  purpose  of  this  thesis  is  to  find  and  investigate  the  opti¬ 
mal  trajectory  which  completes  a  missile  intercept  such  that  some  mea¬ 
sure  of  information  is  maximized.  In  this  way  the  filter  can  make  the 
most  accurate  state  estimate  possible.  It  is  not  suggested  here  that 
this  be  the  only  measure  of  performance,  but  that  it  be  an  element  in 
developing  a  composite  performance  criteria. 

The  particular  missile  under  consideration  is  a  highly-maneuver 
able,  short-range,  bank-to-turn,  air-to-air  missile  capable  of  100  nor 
mal  g's.  The  missile  model,  that  is,  the  aerodynamics,  the  autopilot, 
the  guidance  law,  etc.,  is  completely  developed  in  Reference  2.  The 
geometry  of  the  intercept  and  the  missile  model  are  discussed  in  Sec¬ 
tions  II  and  III. 

The  measure  of  information  chosen  is  the  information  matrix. 
Measures  on  the  information  matrix  imply  the  relative  observability  in 
estimating  the  state  due  to  measurement  functions  alone.  This  means 
that  the  relative  observability  of  the  system  is  obtained  only  from 
the  deterministic  system  model  of  the  dynamics  and  measurement.  The 
information  matrix  and  its  use  as  a  performance  index  is  developed  in 
Section  IV. 

In  Section  V,  the  optimal  control  problem  is  formulated.  The 
parameter  optimization  method  is  then  discussed,  and  all  the  con¬ 
straints  are  derived. 

The  optimization  method  used  to  solve  this  problem  is  the  pen¬ 
alty  function  -  Lagrange  multiplier  method  discussed  by  Fletcher  in 
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Reference  3.  This  is  a  state-of-the-art  method  for  solving  parameter 
optimization  problems.  It  is  particularly  convenient  since  it  requires 
very  little  set  up  time  and  can  be  used  with  equality  constraints  and 
inequality  constraints  at  the  same  time.  Section  VI  contains  a  brief 
discussion  of  this  method  along  with  discussions  of  the  function  eval¬ 
uations  and  numerical  derivatives. 

Finally,  the  results  are  presented  in  Section  VII.  The  various 
launch  scenarios  are  described,  and  some  concluding  remarks  are  pre¬ 
sented  . 


SECTION  II 


GEOMETRY 


This  section  covers  the  problem  set-up.  The  various  reference 
systems  and  their  relationships  are  discussed,  and  the  launch  and  in¬ 
tercept  geometries  are  described. 

2 . 1  Inertial  Axis  System. 

The  inertial  axis  system  is  a  standard  right-handed  system  with 
the  positive  Z^.  axis  pointing  toward  the  center  of  the  earth.  The 
origin  is  at  mean  sea  level  with  the  missile  launch  point  lying  on  the 
negative  Z  axis.  Because  of  this,  the  magnitude  of  the  component 
of  the  respective  position  vector  corresponds  to  the  altitude  of  the 
missile  or  target. 

2.2  Missile  Axis  Systems. 

There  are  two  axis  systems  associated  with  the  missile  -  a  wind 
axis  and  a  body  axis  system.  These  two  systems  share  a  common  origin, 
assumed  to  be  at  the  missile  center  of  gravity. 

The  wind  axis  system  is  a  standard  right-handed  system  with  the 
Xy  axis  tangent  to  the  flight  path  and  positive  forward.  The  Z^  axis 
is  perpendicular  to  the  axis,  in  the  plane  of  symmetry,  and  positive 
down  if  the  missile  is  in  normal,  level  flight.  The  Yy  axis  is  situ¬ 
ated  perpendicular  to  the  X^Z^  plane  to  complete  the  right-hand  sys¬ 
tem.  The.  angular  relationship  with  the  inertial  system  is  through 
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the  velocity  yaw  angle,  x»  the  velocity  pitch  angle,  y,  and  the  veloc¬ 
ity  roll  angle,  y.  This  angular  relationship  is  derived  in  Reference 
4  and  given  by 
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The  missile  body  axis  system  is  described  by  the  positive  Xg 
axis  directed  through  the  nose  of  the  missile.  The  Z_.  axis  is  perpen- 
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dicular  to  the  Xg  axis,  in  the  plane  of  symmetry,  and  positive  down  if 

the  missile  is  in  normal,  level  flight.  The  Y_  axis  is  situated  to 

o 

complete  the  right-hand  system.  It  is  assumed  that  there  is  no  side¬ 
slip  angle;  therefore,  the  body  axes  and  the  wind  axes  are  related 
only  through  the  angle  of  attack, a  .  This  relationship  is  also  derived 
in  Reference  4  and  is  given  by 
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The  missile  body  axis  system  is  related  to  the  inertial  system 
through  the  body  yaw  angle  ip,  and  the  body  pitch  angle  0.  These  two 
angles  completely  define  the  orientation  of  the  body  longitudinal  axis 
in  inertial  space,  and  are  used  as  the  control  variables  in  the  opti¬ 
mal  control  problem.  The  inertial  to  body  angular  relationship  is  de¬ 
rived  by  successive  rotations  about  the  yaw  axis  and  the  pitch  axis 
and  is  given  by 
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The  relationship  of  the  missile  axis  system  to  the  inertial  axis  sys¬ 
tem  is  shown  in  Figure  1. 

The  angles  a  and  y  can  be  found  from  the  angles  ip,  0,  x»  and 
y.  ip  and  0  are  the  control  angles,  and,  therefore,  known;  x  and  Y  are 
part  of  the  state  vector,  and,  therefore,  also  known.  The  solution  is 
found  by  transforming  from  the  inertial  system  through  the  wind  axes 
to  the  body  axes.  If  Equations  (1)  and  (2)  are  rewritten  as 
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the  relationship  between  the  body  axes  and  the  inertial  axes  is  given 
by 
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Performing  the  multiplication  leads  to 
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Equations  (3)  and  (5)  must  be  equal  transformations.  By  setting  the 
components  equal,  the  equations  for  y  and  a  are  found  to  be 
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tanu  =.  cosesin(^-x) _ 

sin0cosy-cos0cos(^-x) siny 


(8) 


tana  =  ^s^nQcos">'~cosQcos^~^)  siny)  cosy+cos0sin(^-x)  siny  , 
co  s  0c  o  s  (  iJj-x)  c  o  s  Y+s  in  0  s  iny 

2.3  Target  System. 

For  this  problem,  the  target  is  simply  assumed  to  be  a  point 
mass  with  position  and  velocity  components  in  inertial  space. 

2.4  Launch  Geometry. 

Figure  2  shows  a  view  of  the  launch  geometry.  At  launch,  the 
missile  is  assumed  to  be  on  the  negative  axis  with  the  body  axis 
system  parallel  to  the  inertial  axis  system.  The  target  initial  posi 
tion  is  determined  by  the  initial  range  and  the  off-boresight  angle, 
BAY, from  the  missile.  Target  direction  is  determined  by  the  initial 
aspect  angle,  AAY.  This  angle  is  defined  to  be  the  angle  between  the 
target  velocity  vector  and  the  line  of  sight  between  the  missile  and 
target.  The  missile  and  target  are  assumed  to  be  co-altitude  and  co¬ 
velocity  at  launch. 

After  launch  the  inertial  axis  system  remains  stationary.  The 
target  remains  at  the  launch  altitude,  but  the  missile  is  free  to  ma¬ 
neuver  in  any  direction.  Figure  3  shows  a  representation  of  the  in¬ 
tercept  geometry.  The  relative  distances  between  the  missile  and  tar 
get  are  also  shown.  The  two  inertial  angular  measurements  are  the 
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azimuth  az  and  the  elevation  el  of  the  target  with  respect  to  the  mis¬ 
sile.  These  angles  are  defined  by  the  relative  distances  as  follows 


az  =»  tan_1(YR/XR) 

el  -  tan-1(-ZR/(X^  +  yJ)5*)  . 


(10) 


SECTION  III 


MODELS 


This  section  defines  the  missile  model  and  all  of  its  compo¬ 
nents,  the  target  model,  and  the  standard  atmosphere. 

3 . 1  Missile  Model. 

The  missile  model  used  is  the  three  degree  of  freedom  version 
of  the  model  derived  in  Reference  2.  Here,  the  missile  is  assumed  to 
be  able  to  instantly  assume  whatever  body  angles  4*  and  0  are  commanded 
by  the  parameter  optimization  routine.  Since  it  is  the  trajectory 
that  is  being  optimized,  it  is  assumed  that  this  trajectory  can  be 
commanded  perfectly  by  an  autopilot;  therefore  there  is  no  autopilot. 
Contributions  made  by  the  control  surface  deflections,  by  sideslip, 
and  by  moments  of  inertia  are  assumed  to  be  negligible. 

The  missile  launch  takes  place  at  t  =  0  in  accordance  with  the 
geometry  in  Section  2.4.  The  missile  is  not  allowed  to  guide  for  the 
first  .4  second,  so  for  this  time  period  the  control  angles  are  zero. 
At  the  end  of  this  delay  the  missile  is  allowed  to  assume  any  position 
commanded  by  the  optimization  routine. 

The  differential  equations  of  motion  for  the  missile  are  as  de¬ 
veloped  in  Reference  4  for  flight  over  a  flat  earth: 


L3 
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x1  =  =  VmcosYcosx 

x2  =  Ym  =  V^cosYsinX 

*3  =  hi  =  'VMsinY 

i.,  =  V  =  (Tcosot  -  D  -  WsinY)/m 

x,.  =  Y  =  ((Tslna  +  L)cosy  -  WcosYJ/V^m 

x^  =  X  =  (Tslna  +  L) siny/ (V^mcosY) 

where 


(ID 


T  »  f(t) 

W  =  f(t) 
m  =  W/g 

D  *  fCa»VZM) 

L  -  f(a,VM,ZM) 

a  -  f(X.Y.'M) 
y  -  f(X.Y.'M) 

g  =  acceleration  due  to  gravity,  assumed  constant. 


(12) 


The  functions  for  a  and  y  were  discussed  in  Section  2.  The  models  for 
thrust,  weight,  drag  and  lift  remain  to  be  developed. 

3.1.1  Thrust  and  Weight  Models.  The  missile  weighs  165  lb  at  launch. 
The  propellent  (50  lb)  burns  at  an  assumed  constant  mass  flow  rate  for 
2.6  seconds.  During  the  burn,  the  thrust  is  assumed  to  be  constant, 
4711.5  lb.  There  is  no  other  contribution  to  weight  loss  during  the 
missile's  time  of  flight.  The  thrust  and  weight  profiles  are 
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described  by  the  expressions 


T  = 


W  = 


4711.5  lb 
0 

(165  -  19.23t)  lb 
115  lb 


0  <  t  <  2.6  sec 

2 . 6<  t 

0  <  t  <  2.6  sec 

2 . 6<  t 


(13) 

(14) 


3.1.2  Aerodynamics  Model.  The  components  of  the  missile  aerodynamic 
forces  along  the  and  axes  are  given  by  the  equations 


Fy  =  -(C  +  Cv  a  +  C  .a2)  \  p  V2  S 

X  a  X  X  Z  ^  n 

o  a  a 


FZ--<C»n  +CNa)I0VMS 
oT  a 


(15) 


where  S  is  the  missile  reference  area.  The  atmospheric  density  p  is 
calculated  using  a  standard  atmosphere  model. 

The  aerodynamic  coefficients  in  Equations  (15)  are  found  by 
table  look-up  and  linear  interpolation.  The  values  are  listed  in 
Tables  1  and  2.  These  tables  and  the  interpolation  routines  have  been 
taken  directly  from  Reference  2.  Drag  and  lift  are  then  calculated 
from  Equations  (15)  in  accordance  with  Figure  4. 

D  *  cosa  -  Fz  sina 

L  -  -F_  cosa  +  Fv  sina 


(16) 


TABLE  2.  AERODYNAMIC  COEFFICIENT,  C  (ct,M) 


Figure  5.  Boresight  Angle 
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3.1.3  Seeker  Model.  The  actual  seeker  and  how  it  works  is  not  model¬ 
ed.  However,  to  simulate  the  presence  of  a  seeker,  a  limit  is  imposed 
on  the  angle  between  the  positive  Xg  axis  and  the  line  of  sight  to  the 
target.  This  angle,  shown  in  Figure  5,  is  called  the  boresight  angle 
and  is  limited  to 

|n|  <  60°.  (17) 

This  limit  represents  the  maximum  angle  off  boresight  that  the  seeker 
can  track  the  target.  It  is  a  physical  limitation;  if  the  boresight 
angle  were  to  become  greater  than  60°,  the  seeker  would  break  lock. 

The  boresight  angle  is  found  by  using  the  definition  of  a  dot 
product.  From  Figure  5,  it  is  seen  that 

|0A|  | AT |  cosft=0A  ’  AT  (18) 

where  0  represents  the  missile  center  of  gravity,  and  A  is  the  seeker 
gimbal  pivot  point.  The  inertial  axis  system  shown  is  actually  a  par¬ 
allel  system  through  the  missile  c.g.  T  represents  the  target  posi- 
ion. 

A  unit  vector  in  the  direction  ^5a  can  be  defined  using  the 
Euler  angles  from  Equation  (3),  that  is, 

OA/jOAl  =  (ijCos9cos^  +  j^cos9sin^  -  k^sin9) .  (19) 


i 
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The  vector  AT  is  the  relative  range  vector 

"  •  Vi  +  Vi  +  zk  *1  <20) 

where  ^ 

directions  respectively  (see  Figure  3)  .  The  boresight  angle  is  then 
found  to  be 


fR,  and  ZR  are  the  relative  ranges  in  the  X^,  , 


and  Z, 


cosfi  = 


XR  cos9cosi|J  +  Yr  cos0sin^ 


-  ZR  sin8 


<XE2  +  YE 


2.  U 

+  V 


(21) 


3.1.4  Accelerometer  Model.  The  normal  acceleration  is  limited  to  a 
maximum  of  100  g's.  Missile  normal  acceleration  is  modeled  by 

a7  -  F_  /  W  (22) 


where  Fz  is  the  component  of  aerodynamic  force  along  the  ZR  axis  as 
defined  in  Equations  (15),  and  W  is  the  weight  from  Equation  (14). 
The  limit,  for  structural  reasons,  is 


<  100  g's 


(23) 
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3.2  Target  Model. 

The  target  is  considered  to  be  a  non-maneuvering,  non-acceler¬ 
ating,  point  mass  moving  through  inertial  space.  Therefore,  the  tar¬ 
get  model  consists  of  only  the  three  inertial  position  components 

Xt  =  *r  +  vt  c 
0  x 

Y  -  Y  +  VT  t  (24) 

0  y 


For  this  problem,  the  velocity  components  are  as  follows: 

VT  cos(AAY  +  BAY) 

VT  s in ( AAY  +  BAY)  (25) 

0 

where  VT  is  the  constant  target  velocity  and  AAY  and  BAY  are  the  ini¬ 
tial  aspect  angle  and  of f-boresight  angle,  respectively,  as  defined  in 
Section  2.4. 

3.3  Atmospheric  Model. 

The  atmospheric  model  is  needed  to  compute  the  speed  of  sound 
and  the  density  for  a  given  altitude.  This  model  is  a  standard,  con¬ 
stant  gravity  atmosphere  for  the  troposphere  (0  to  36089  ft)  and  the 
stratosphere  (36089  to  82021  ft) .  It  is  the  same  model  developed  in 


Reference  2. 


SECTION  IV 


THE  INFORMATION  MATRIX 


In  Reference  1,  a  dynamical  system  used  in  an  onboard,  optimal 
filter  to  estimate  relative  position,  relative  velocity,  and  target 
acceleration  is  presented.  In  order  to  measure  the  observability  of 
the  state  estimate  which  can  be  made  from  this  dynamical  system  and 
measurement  function,  the  information  matrix  is  used.  In  this  section, 
the  information  matrix  and  measures  associated  with  the  information 
matrix  are  developed  as  the  performance  criteria  for  the  optimization 
problem. 

4.1  Definition. 

The  information  matrix  is  defined  by  Reference  5  as 

I(t,t0)  =  /  $T(x,t)HT(x)V_1(x)H(x)  <Kx,t)dx  (26) 

where  <Hx,t)  is  the  state  transition  matrix  used  by  the  filter.  H(x) 
is  the  matrix  of  partial  derivatives  of  the  measurement  function  with 
respect  to  the  state  variables  evaluated  on  the  nominal  path,  and 
V(x)  is  the  measurement  noise  power  spectral  density.  Each  of  these 
elements  will  be  discussed  in  this  section. 

The  information  matrix  is  related  to  the  inverse  of  the  error 
covariance  matrix  P  ^(t).  If  it  is  assumed  that  there  is  no  process 
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noise, 

by 


for  any  time  period 


ci t0  V 


the  information  matrix  is  given 


1  ( t j ’ t i>  ~  f  ~  $  (t0’tj)  po  <i>(to,tj) 


(27) 


where  P  is  the  a  priori  error  covariance  (Reference  6) .  If  there  is 

0 

no  a  priori  information,  i.e.,  if  PQ^  =  0, 

I(t.,t.)  -  P'1Ct.).  (28) 

J  i  J 

The  larger  the  eigenvalues  of  the  information  matrix,  the  smaller  the 
eigenvalues  of  the  corresponding  error  covariance  matrix  will  be,  and 
the  more  precise  the  estimate  will  b«.  This  is  the  basis  for  using 
the  information  matrix  as  the  performance  index.  The  trajectory 
which  maximizes  the  eigenvalues  of  the  information  matrix  is  also  the 
trajectory  which  gives  the  jpost  precise  estimate. 

A  scalar  is  needed  for  the  performance  index.  The  two  opera¬ 
tions  available  are  the  trace  and  determinant.  The  trace  is  equal  to 
the  sum  of  the  eigenvalues  of  a  matrix;  the  determinant  is  equal  to 
the  product  of  the  eigenvalues.  Because  of  its  simplifying  charac¬ 
teristics,  primarily  that  the  trace  and  integration  operations  are 
interchangeable,  the  trace  of  the  information  matrix  has  been  chosen. 
Since  this  is  an  optimization  problem  solved  on  a  finite  word  length 
computer,  a  constant  weighting  matrix  is  included  to  normalize  the 
terms  in  the  trace.  This  weighting  matrix  does  change  the  eigen¬ 
values,  but  since  it  is  a  constant,  the  maximizing  has  the  same  effect 
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The  performance  index  now  becomes 

J  =  tr  W  I  (tf,tQ) 

=  tr  W  /  f<t  T(x,t  )HT(X(t))V_1(t)H(X(t))  <Kx,t  )dx.  (29) 

co 

Since  the  trace  and  the  integral  operations  are  interchangeable  and 
the  weighting  matrix  is  a  constant,  J  can  be  rewritten  as 

J  =  /  ftr(W$T(x,  tf)HT(X(x))v“1(x)H(X(x))$(x,t  ))dx.  (30) 

t0 

This  defines  the  form  of  the  performance  index  used  for  this  problem. 
The  elements  of  this  equation  remain  to  be  defined. 

4.2  The  Dynamical  System. 

The  dynamical  system  used  by  the  filter  (Reference  1)  is  of  the 

form 


X(t)  =  F  X(t)  +  G  u(t) 

(31) 

N 

rt 

w 

II 

0Q 

/~N 

X 

rt 

V— ✓ 

V-/ 

+ 

< 

s 

rt 

• 

(32) 

The  dynamics  equation  (Equation  31)  is  linear  in  rectangular  coordi¬ 
nates.  X  is  a  nine-state  vector  of  relative  position,  relative  veloc¬ 
ity  and  target  acceleration  given  by 
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X  = 


yr  zr 


03) 


where 


XR  -  XR(t) 

yr  *  Vc) 

ZR  '  ZR(t> 
UR  *  "R(t) 

VR  '  V'> 

»R  ‘ 


-  xT(t)  -  X^t) 

*  V0  ‘  VM(t) 

'  ZT<C>  '  ZM(C) 

‘  \(t>  '  \<£> 

'  \M  -  \<£> 

”  VT  -  VM  <C>> 

z  z 


(34) 


and  a  ,  a  ,  and  a  are  the  inertial  target  acceleration  components. 
X  Y  Z 

Equation  (31)  in  scalar  form  is 


XR  =  UR 

Y  =  V 
R  R 

XR  =  WR 


Ur  =  X  '  X 


VR  aTy  “  \ 


WR=aTz'X 

S '  'x  % 

iT  =  -x  a 

Y  Y 

a  =  -X  a 
Z  Z 


(35) 


he 
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where  A  is  a  target  acceleration  constant.  The  control,  u(t),  is 
assumed  known,  and  consists  of  the  three  missile  inertial  acceleration 


and 


aM„ 


components,  a^,  aHy’ 
the  matrix  F  from  Equation  (31) 


From  Equation  (35)  it  can  be  seen  that 
is 


0 

1  I 

-  0  1 

-  - 

F  = 

0 

1  0 

1  I 

1“  “ 

_  0 

1 

i  0 

I 

i-A  _ 

(36) 


and  the  matrix  G  is 


G 


0  I  0  I  0 
0  I -I  i  0 
0  l  0  I  0 


(37) 


where  each  partition  is  3x3,  I  is  the  identity  matrix,  and  A  is  a 
diagonal  matrix  of  the  target  acceleration  constants. 

The  measurement  equation  (Equation  (32)  )  is  nonlinear  in  rectan¬ 
gular  coordinates.  The  measurement  z  is  a  two-dimensional  angle  mea¬ 
surement,  g(X)  is  the  measurement  function,  and  v(t)  is  a  zero  mean, 
white  noise  process  associated  with  the  measurement,  with  power  spectral 
density  V(t).  The  measurement  function  consists  of  the  two  inertial 
angle  measurements  for  azimuth  and  elevation  (see  Figure  3) 


g(X(t)) 


tan"1 (Yr/Xr) 
tan_1(-ZR/(XR2  +  Yr2)^) 


(38) 
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The  linearized  measurement  equation  is 


<5z(t)  =  H(X(t)  )SX  +  v(t). 


where  H(X(t))  is  the  partial  of  g  with  respect  to  X  evaluated  on  the 
nominal  path.  The  elements  of  H  are 


8ly  S1 

H(X(t) )  =  X1  X2 


0  0  0  0  0  0  0 


g?  g2  g2  0  0  0  0  0  0 

xx  x2  x3 


where 


H11  81„  R 

X^  az 


H12  “  glY  R 

X2  az 


=  ¥r._l 


“21  '  82x,  ‘  R2,  Raz 

1  el 

yrzr  .  _1_ 

22  §2X  R2  Raz 

X2  R  el  az 


l23  =  g2. 
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where 

Saz  '  *  <XR2  +  ^ 

*.l  *  Eel(t)  '  «R2  +  XR2  +  • 

as  shown  in  Figure  3. 


(42) 


4.3  The  State  Transition  Matrix. 

Because  of  the  simplicity  of  the  matrix  F  (Equation  (36)  ) , 
the  state  transition  matrix  has  been  derived  in  closed  form.  The 
elements  of  3>(t,r)  are  given  by  (Reference  1) 


$,  .  0  0  ft .  0  0  $..,  0  0 

11  14  17 

0  0  0  $25  0  0  $>g  0 

0  0  $33  0  0  $36  0  0  *39 

o  o  o  $44  o  o  0  0 

$(t,T)  =  0  0  0  0  $55  0  0  $58  0 

0  ooo  0  $66  0  o  $69 

0  0  0  0  0  0  $77  0  0 

0  0  0  0  ooo  ^80 


(43) 


0 


ooo 


ooo 


where 


$n(t,T)  =  $22(t,t)  =  4>33(t,i)  =  1 

$44(t,x)  »  $55(t,T)  =  $66(t,T)  -  1 

$14(t,r)  *  $25(t,r)  =  $36(t,x)  *  t  -x  =  At 

<J>77(t,x)  *  <J>g8(t,x)  =  $9g(t,x)  =  exp(-X  At) 

4>47(t,x)  =  $58(t,x)  =  $69(t,x)  =  -  y(exp(-X  At)  -1) 

<5*17 <t: ,T)  =  «28(t,x)  =  $39(t,x)-  =  -ij(exp(-X  At)  +  XAt  -  1) 

X 


(44) 


and  where  X  is  the  target  acceleration  constant  (Reference  1). 

Since  the  information  matrix  requires  $(x,t),  the  identity  for 
the  inverse  of  a  state  transition  matrix 

<Kx,t)  =  $-1(t,x)  (45) 


can  be  used.  Hence,  the  elements  of  <J>(x,t)  are 


^ll(x,t)  *  $22^T,t^ 

$44(x,t)  =  $55(x,i) 
^(x.t)  *  $25(x,t) 

$77(x,t)  *  $88(t,t) 

<j>47  (x  «t)  =  $3g(x,t) 
<J>17(x,t)  =  ^(x.t) 


<J>33(x,t)  =  1 

$66^T,t^  =  1 

$36(x»t)  =  x  -  t  »  -At 

<f>gg(x»t)  =  exp(AAt) 

$6g(x,t)  =  -y(exp(x  At)  -  1) 

1 

$3g(x,t)  ■  •—( exp ( \  At)  -X  At  -  1) 


(46) 
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4.4  Angle  Measurement  Noise. 

The  matrix  V(t)  is  a  2x2  matrix  representing  the  power  spectral 
density  of  the  angle  measurements  due  to  certain  target-sensor  charac¬ 
teristics,  that  is. 


V(t) 


vu(t>  0 

0  V22(t) 


VU(t) 


+  b 


az 


V22^C)  2  +  b  • 
el 


(47) 


V^^(t)  denotes  the  power  spectral  density  associated  with  measuring 
azimuth  and  V ^ (t)  denotes  the  power  spectral  density  associated  with 
measuring  elevation.  The  term  a  is  a  constant  associated  with  the 
received  signal  power  noise.  The  term  b  represents  the  uncertainty 
in  target  position,  possibly  due  to  atmospheric  refraction  of  the 
signal  to  the  sensor  and  fading.  For  this  problem, 

2  2 

a  =  .25  rad  ft  sec 

-8  2  (48) 
b  ■  56.25  x  10  rad  sec 


are  typical  values  for  current  sensors. 
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V(t)  may  be  simplified  by  factoring  out  a,  that  is. 


V(t)  =  a 


0 


(49) 


where  c  =»  b/a. 


4 . 5  Weighting  Matrix. 

The  weighting  matrix  W  is  a  constant,  9x9  diagonal  matrix 
used  to  give  relative,  equal  values  to  the  terms  in  the  information 
matrix.  The  off-diagonal  terms  are  all  zero,  while  the  diagonal  terms 
are  defined  as  follows: 


Uu  *  “22  '  W33  ■  W1  ■  1 
“44  *  “55  *  "66  '  W2  "  Cf 


-2 

f 

-2t, 


W77  *  W88  *  W99  ’  W3  “  e  f* 


(50) 


Wl’  W2’  an<*  ^3  weiSht  t^ie  positi-00  terms,  the  velocity  terms  and  the 
target  acceleration  terms,  respectively.  They  have  been  chosen  in  an 
effort  to  get  the  maximum  numerical  value  of  each  of  these  terms  to  be 
approximately  equal  to  one.  Without  this  relative  weighting,  the  tar¬ 
get  acceleration  terms  completely  dominate  the  other  factors.  This 
becomes  obvious  in  Section  4.6  where  the  performance  index  is 
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simplified,  and  the  effects  of  the  weights  can  be  seen. 

4.6  The  Performance  Index. 

After  multiplying  out  the  terms  in  Equation  (30)  and  performing 
the  trace  operation  inside  the  integral,  J  becomes 


J  -  /Cf((4u  HH1«11+  HH2  *22  +  *33  HH3  *33)W1 

co 


14 

HH1 

*14  + 

$ 

w25 

Hn2 

$  + 
z5 

4 

36 

HH3 

*36)W2 
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where  $  =  <J>(T,tf),  and  where 


HH1  =  HH1(t) 


hh2  -  hh2(t) 


hh3  =  hh3(t) 
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X3 
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(52) 


The  performance  index  further  simplifies  to 


SECTION  V 


PROBLEM  FORMULATION 


This  section  formulates  the  problem  as  an  optimal  control  pro¬ 
blem  and  presents  all  of  the  constraints-  The  problem  is  then  con¬ 
verted  to  a  parameter  optimization  problem. 

5.1  The  Optimal  Control  Problem. 

The  objective  is  to  find  the  control  which  maximizes  the  infor¬ 
mation,  completes  the  air-to-air  intercept  (zero  miss  distance) ,  and 
satisfies  all  of  the  dynamical,  physical,  and  aerodynamic  constraints, 
that  is. 


Maximize  Equation  (53)  (54) 

subject  to  the  terminal  constraint  of  zero  miss  distance  and  the  sys¬ 
tem  model  developed  in  Section  3.  The  dynamical  part  of  this  model 
takes  the  form  of  differential  constraints  and  is  nonlinear  in  the  form 

x  =  f  (t ,  x,  u)  (55) 

where  the  control  u  consists  of  the  body  angles  and  0. 

The  initial  conditions  x^  for  the  differential  equations  of  the 


missile  are 
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xx(0)  =  0 
*2(0)  =  0 
*3(0)  -  -h 

V°>  *  VM 
x5(0)  =  o 

x6(0)  =  0 


0 

0 


(56) 


where  the  values  h„  and  V__  are  input  values  for  initial  altitude  and 
0  M0 

velocity. 

The  terminal  constraint  of  zero  miss  distance  is  formulated  as 
three  equality  constraints  in  the  form 


“V 

r(xT  -  >  1  h:  i 

Tf  nf  Af 

v  (tf,Xf)  = 

^2 

= 

(Y  -  Y  )  /  Y 
'  T  M  J  T 

f  f  f 

*3 

(Z  -  Z  )  /  z 
v  Tf  nf '  '  Tf 

(57) 


where  the  values  of  ¥  denote  the  actual  constraint  residuals  for  the 
current  control.  The  object  is  to  drive  these  values  to  zero.  X^.,  YT, 
and  Z,p  are  the  target  inertial  position  components.  The  subscript  f 
denotes  the  condition  at  the  final  time. 

The  remaining  constraints  are  due  to  physical  or  aerodynamic 
limitations  and  are  formulated  as  inequality  constraints,  0(t ,x,u)>  0. 
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The  limitations  are 

t^  <  8  sec 

| a  |<  100  g's 
ZB  ‘ 

|Q|  <  60°  (58) 

|a|  <  20° 

|y|  <  80° 


The  final  time  is  limited  to  no  more  than  8  seconds  in  order  to  ensure 
that  the  missile  still  has  enough  velocity  to  maneuver.  The  missile 
normal  acceleration  is  limited  to  100  g's  for  structural  reasons.  The 
missile  seeker  has  a  maximum  angle  limit  of  60°  off  boresight,  hence 
the  constraint  on  boresight  angle.  Maximum  angle  of  attack  is  20°,  as 
limited  by  the  aerodynamics  tables.  The  limit  on  the  flight  path 
angle  is  included  to  keep  the  equations  of  motion  defined.  The  value 
of  x  (Equations  (11)  )  becomes  infinite  at  y  =  90°.  Instead  of  re¬ 
defining  the  equations,  a  state  inequality  constraint  is  used  to  keep 
y  away  from  90°  because,  it  was  felt,  the  final  trajectory  would  not 
require  Y  above  80°.  Until  the  converged  solution  is  found,  these  in¬ 
equality  constraints  may  or  may  not  be  satisfied. 

During  the  optimization  process,  the  inequality  constraints  are 
evalutated  according  to  their  constraint  residuals.  For  the  final 
time  constraint 


1  -  tf  /  8 


(59) 
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The  remaining  inequality  constraints  could  be  violated  at  any  time 
from  tg  to  t^.  Therefore,  the  constraint  residuals  are  evaluated  as 
integral  constraints,  that  is. 


The  performance  index  (Equation  (54)  )  and  the  integral  con¬ 
straints  are  formulated  as  differential  equations  with  zero  initial 
values  and  evaluated  along  with  the  state  differential  equations. 

This  augments  the  state  vector  (Equation  (11)  )  by  five:  x^  for  the 
performance  index,  and  Xg  through  x^-^  for  the  integral  constraint 
residuals.  The  values  of  the  performance  index  and  the  integral  con¬ 
straint  residuals  (Equation  (60)  )  are  the  values  of  the  augmented 
state  variables  at  the  final  time,  that  is, 


Q 


3 


x 

9f 


(61) 
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5.2  The  Parameter  Optimization  Method. 

The  preceding  forms  an  optimal  control  problem  to  find  the  con¬ 
trol  u(t)  which  maximizes  the  performance  index.  To  solve  this,  the 
problem  is  converted  to  a  parameter  optimization  problem  where  the 
constraints  are  handled  by  penalty  functions. 

The  problem  of  having  a  free  final  time  is  converted  to  a  fixed, 
normalized  final  time  by  defining  a  new  independent  variable. 

T  -  t  /  tf  (62) 


or 


dt  *  tf  dT. 

Thus,  the  equations  of  motion  become 


dx  =  dx 
dt  “  tfdx 


f (tfT,x,u) 


(63) 


(64) 


or 


dx 

dx 


g(x,x,u,tf)  , 


(65) 


where  x^  =  0,  x^  =  1.  Final  time  is  now  a  parameter  to  be  found  along 


with  the  control.  The  control  is  now  a  function  of  X,  the  new 
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Independent  variable. 

The  continuous  control  function  is  approximated  by  two  cubic 
splines,  one  for  each  control,  with  unknown  parameters,  that  is, 

u  =  u(a,x)  ,  (66) 

where  a  is  a  vector  of  n  unknown  parameters 

»„].T  (67) 

For  each  control  function  these  parameters  are  the  values  of  the  con¬ 
trol  at  T  =  0,  T  =  h,  t  =  1,  and  the  slope  at  each  end  point.  For  any 
value  of  a,  a  cubic  spline  is  fit  to  the  parameters,  and  the  control 
can  be  found  for  any  T  between  0  and  1. 

The  parameters  to  optimize  now  are  the  final  time  and  the  para¬ 
meters  of  the  cubic  splines  which  define  the  control  function.  They 
can  be  lumped  together  as  one  vector  of  length  n  +  1 

b  =  Gtf  aT*  (68) 

From  Equations  (65),  (66),  and  (68),  it  can  be  seen  that  x^  =  x^(b) . 
Therefore,  since  the  performance  index  is  a  function  of  x^  and  t^ ,  J 
can  now  be  written 


J  -  0(b) 


(69) 


where  4>( b)  =  -x 


The  problem  can  now  be  restated  as 


Minimize  J  =  <p(b) 


subject  to 


nb)  =  0 

0(b)  >  0  , 


where 


‘■Qf  »]T 

u  =  u(a,x) . 
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SECTION  VI 


SOLUTION  METHOD 


This  section  gives  a  brief  description  of  the  problem  solution 

method . 

6.1  Input . 

The  input  required  for  this  program  is  initial  missile  and  tar¬ 
get  positions  and  an  initial  set  of  parameters,  b.  The  missile  and 
target  are  set  up  as  described  in  Section  2.  An  initial  altitude, 
velocity,  range,  of f-boresight  angle  and  aspect  angle  must  be  speci¬ 
fied.  The  control  is  represented  by  two  cubic  splines,  one  for  each 
control  angle.  Each  spline  function  is  assumed  to  have  5  parameters  - 
the  control  angles  at  T  =  0,  T  =  *5,  and  at  t  =  1,  and  the  slope  at 
each  end  point.  Since  t^  is  also  a  parameter,  a  total  of  11  parame¬ 
ters  must  be  input. 

6.2  Optimization  Method. 

The  optimization  method  used  is  a  penalty  function-Lagrange 
multiplier  method  which  is  completely  developed  in  Reference  3.  The 
penalty  function  used  is 

P(b,B,S)  =  <Kb)  +j  r(b,g)T  S  r(b,B).  (73) 

Here,  b  is  the  n  +  1  vector  from  Equation  (69)  which  contains  the 
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parameters  to  be  optimized.  T(b)  is  an  m  vector  where  m  is  the  total 
number  of  constraints,  that  is, 


ro>,B) 


' 

C1(b)  -Bi 
min(Ci(b)  -&±t  0) 


,  1  5  i  5  k 

,  k+1  1  i  -  m 


(74) 


where 


C 


1 


C(b)  = 


(75) 


C  =  0 
m  m-k 


Hence,  the  first  k  elements  of  C(b)  are  the  equality  constraint  resid¬ 
uals,  and  elements  k+1  to  m  are  the  inequality  constraint  residuals. 

S  is  an  m  x  m  diagonal  weighting  matrix  with  elements  a^,  where  i  goes 
from  1  to  m.  The  m  vector  3  is  used  to  allow  convergence  without  the 
necessity  of  forcing  S  to  infinity. 

This  is  an  iterative  method  where  each  iteration  involves  mini¬ 
mizing  P(b,8,S)  for  a  fixed  8  and  S  using  a  variable-metric  method  for 
unconstrained  optimization.  After  each  iteration,  0  and  S  are  varied 
in  such  a  way  that  the  parameters  tend  to  the  constrained  solution. 

The  aim  is  to  vary  8  so  that,  for  a  constant  S,  the  value  of  the  product 
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0  S  tends  to  the  vector  of  Lagrange  multipliers  for  the  problem.  The 
values  are  only  increased  in  case  the  rate  of  convergence  of  the 
corresponding  C^(b)  to  zero  is  not  sufficiently  rapid.  Convergence  is 
obtained  when  |c^(b)jie^,  for  l^ilm. 

The  computer  subroutine  incorporating  this  method  was  taken 
from  the  Harwell  library  and  is  named  VF01A.  To  use  VF01A,  the  user 
must  supply  a  routine  called  VF01B  which  evaluates  4>(b) ,  Vcf) (b) ,  C(b)  , 
and  VC(b) ,  where 

V<Kb)  =  (76) 

3C. (b) 

VCij(b)  =  3b^-  (77> 


for  1  i  i  -  (n+1)  and  1  1  j  1  m. 

6.3  Function  Evaluation. 

The  function  evaluation  consists  of  finding  the  performance 
index  J  =  4>(b)  ,  and  the  constraint  violations,  C(b)  .  These  values 
are  then  used  to  form  Equation  (73)  .  Figure  6  is  a  block  diagram 
showing  the  elements  of  the  program  which  handle  the  function  evalua¬ 
tion.  Each  block  represents  a  subroutine. 

Subroutine  SG  is  called  from  VF01B  with  the  current  set  of 


parameters  b.  SG  is  the  driving  routine  which  sets  all  the  initial 
conditions,  calls  the  integrator,  and  then  returns  the  values  of  4> 
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and  C.  SPLINE  is  called  Immediately  upon  entry  to  SG  to  evaluate  the 
coefficients  of  the  cubic  spline  functions  of  the  control  for  a  fixed 
b.  The  control  can  then  be  evalutated  explicitly  as  a  function  of  T. 

The  integrator  RUNGE  is  a  fixed-step,  fourth-order  Runge-Kutta 
method.  The  system  differential  equations,  the  performance  index, 
and  the  integral  constraints  are  all  integrated  simultaneously. 

The  complete  system  model  is  contained  in,  or  called  by  DERIV. 
For  every  time  t,  DERIV  finds  the  control  by  calling  SPLINE  and  evalu¬ 
ates  the  differential  equations.  There  are  11  differential  equations: 
6  for  the  state  equations,  1  for  the  performance  index,  and  4  for  the 
integral  constraints.  The  aerodynamics  tables  required  to  evaluate 
drag  and  lift  are  contained  in  AERO.  ATMOS,  the  standard  atmosphere 
model,  must  be  called  to  find  the  Mach  number  and  the  atmospheric  den¬ 
sity  for  the  instantaneous  altitude  and  velocity.  The  aerodynamic 
tables  (Tables  1  and  2)  are  functions  of  Mach  or  of  Mach  and  angle  of 
attack.  The  routines  used  to  linearly  interpolate  the  aerodynamic  co¬ 
efficient  tables  are  INTRP1  and  INTRP2,  respectively. 

6.4  Numerical  Derivatives. 

To  obtain  V$(b)  and  VC(b)  as  required  for  VF01A,  numerical  de¬ 
rivatives  are  evaluated  by  the  method  of  central  differences.  For  any 
function  G(x)  ,  the  numerical  derivative  with  respect  to  the  vari¬ 
able,  x^»  is  found  by  the  formula 

3G  _  G(x^  +  ~  *^xk  ” ^  xk^ 

9x, 
k 


(78) 
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where 


(79) 


atimization  Algorithm. 


The  algorithm  used  for  constrained,  numerical  optimization  in 


VF01A  is  stated  as  follows: 


1.  Guess  b. 

2.  Set  initial  values  of  S,|3. 

3.  Minimize  P(b,f3,S)  with  respect  to  b. 

a.  Obtain  <J>  and  C  by  calling  function  evaluation  routine. 

b.  Obtain  V((>  and  VC  by  calling  numerical  derivative 
routine. 

c.  Evaluate  P(b,8,S)  and  3P (b, B,S) /9b. 

d.  Use  variable  metric  method  for  unconstrained  minimiza¬ 
tion. 

4.  If  the  final  conditions  are  set  to  desired  accuracy,  go 
to  6. 

5.  Vary  8,  and  S  if  necessary;  to  to  2. 

6.  Return  to  executive  program. 


SECTION  VII 


RESULTS 

In  this  section  the  launch  scenarios  are  described,  and  the 
results  are  discussed. 

7.1  Launch  Scenarios. 

Three  different  launch  scenarios  have  been  chosen  for  the  com¬ 
putation  of  optimal  information  trajectories.  Each  intercept  is  ini¬ 
tiated  at  10,000  ft  altitude  with  the  missile  and  target  co-altitude 
and  co-speed  at  .9  Mach.  The  scenarios  differ  in  the  launch  range 
and  aspect  angle  only,  as  the  boresight  angle  is  zero  for  all  engage¬ 
ments.  The  engagements  are  initiated  at  3000  ft,  0°  aspect;  3000  ft, 
30°  aspect;  and  at  7000  ft,  60°  aspect.  The  target  for  each  intercept 
is  flying  straight  and  level  and  nonaccelerating. 

In  order  to  test  the  system  model  and  for  comparison  purposes, 
one  additional  intercept  has  been  run,  that  of  minimum  final  time. 

This  is  done  simply  by  replacing  the  performance  index  with 

J  “  t^  .  (80) 

Everything  else  remains  the  same.  The  launch  for  this  intercept 
occurs  at  3000  ft  and  30°  aspect  angle. 
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7.2  Results  and  Conclusions. 

Since  the  measurements  for  the  dynamical  system  are  angle  mea¬ 
surements  only,  it  is  expected  that  to  increase  the  information  the 
missile  needs  to  fly  a  path  which  keeps  these  angles  changing.  The 
optimal  trajectories  found  with  this  method  are  consistent  with  this 
expectation.  During  the  intercept,  the  missile  "fishtails"  and  "por¬ 
poises"  to  the  maximum  extent  possible  without  violating  any  con¬ 
straints.  The  angle  rates  for  the  boresight  angle  and  for  the  mea¬ 
surement  angles  are  kept  high  during  the  entire  intercept.  The  bore- 
sight  angle  tries  to  go  from  limit  to  limit. 

Figures  8  through  24  present  the  results  in  graphical  form  for 
each  of  the  intercept  scenarios.  The  missile  launch  point  is  where 
Xj  =  0,  Yj  *  0,  and  Zj  =  -10,000  ft.  Since  nothing  can  happen  for  the 
first  .4  seconds,  the  plots  show  only  the  time  period  for  .4  <t  <  t^. 
The  figures  are  grouped  by  intercept,  each  group  containing  all  the 
results  for  each  scenario.  The  first  two  figures  in  each  group  show 
the  actual  intercept  and  target  trajectories  projected  onto  the  X^-Y^ 
and  the  X^-Z^  planes  respectively.  For  the  horizontal  trajectories, 
the  viewer  is  above  the  X^-Y-j.  plane  looking  down;  for  the  vertical 
trajectories,  the  viewer  is  on  the  positive  axis  side  of  the  X^-Z^ 
plane.  The  next  four  figures  are  for  the  boresight  angle,  the  bore- 
sight  angle  rate,  the  azimuth  and  elevation  (measurement)  angles,  and 
the  azimuth  and  elevation  angle  rates,  all  plotted  versus  time. 

The  next  6  figures  (Figures  25  through  30)  present  the  same  in¬ 
formation  for  the  minimum  time  trajectory.  These  are  shown  primarily 
for  comparison  with  the  best  information  trajectories,  but  they  also 


demonstrate  the  validity  of  the  system  model  and  the  optimization  pro¬ 


cess.  The  minimum  time  results  are  just  what  you  would  expect.  The 
trajectory  is  an  easy,  but  direct  turn  onto  an  intercept  path  which 
approaches  a  straight  line.  There  is  no  variation  in  elevation.  An 
abrupt,  hard  turn  would  induce  more  drag  and,  thus,  make  a  longer  in¬ 
tercept.  For  contrast,  the  best  information  trajectory  is  also  shown 
on  the  minimum  time  trajectory  plots. 

Table  3  lists  the  solved  intercept  problems  by  engagements. 

Each  problem  has  converged  satisfactorily,  and  the  miss  distances  are 
all  within  1  foot  of  the  point  mass  target.  The  only  inequality  con¬ 
straint  which  acts  on  the  converged  solution  is  the  maximum  boresight 
angle.  The  values  of  the  performance  index  are  also  listed,  but  the 
only  ones  that  can  be  compared  are.  the  ones  with  the  same  initial  geo¬ 
metry.  For  the  3000  ft,  30°  launch,  the  performance  index  for  the 
maximum  information  trajectory  is  much  higher  than  for  the  minimum 
time  trajectory,  as  it  should  be. 

Note  the  high  number  of  function  evaluations  and  the  amount  of 
computer  time  required  for  the  best  information  solutions.  Because  of 
this,  it  would  be  impossible  to  use  this  optimization  process  in  an 
onboard  computer  for  a  guidance  scheme.  If  something  like  this  were 
desirable,  however,  it  may  be  possible  to  streamline  the  process  in 
some  way.  This  is  a  possible  area  for  further  research.  For  example, 
the  aerodynamics  tables  could  be  replaced  by  some  approximating  func¬ 
tion,  eliminating  the  table  look-up  and  interpolation  process. 

Another  possibility  is  to  formulate  an  empirical  formula  for 
the  guidance  law  which  would  approximate  the  best  information 


i 
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trajectory.  The  plots  for  the  three  intercepts  show  a  definite  simi¬ 
larity  between  the  scenarios.  On  each  intercept  the  missile  appears 
to  be  attempting  to  do  the  same  thing,  limited  only  by  its  requirement 
to  complete  the  intercept.  If,  after  running  more  trajectories,  this 
is  found  to  be  the  case,  then  an  empirical  formula  might  be  found. 

Another  possible  area  for  future  research  on  this  project  is  to 
find  out  how  the  trajectories  would  differ  if  some  information  about 
the  state  is  removed.  For  example,  if  in  Equation  (53)  is  set  to 
zero,  the  best  information  trajectory  without  target  acceleration  in¬ 
formation  would  be  found.  This  is  because  the  target  acceleration 
terms  in  the  state  transition  matrix  would  be  blanked  out.  If,  it  is 
found  that  the  optimal  trajectory  is  satisfactory  without  the  target 
acceleration  information, then  this  information  would  not  have  to  be 
estimated  by  the  onboard  filter  for  the  guidance  law. 


Figure  15.  Best  Information  Boresight  Angle 
3000  Feet,  30  Degrees. 


Figure  16.  Best  Information  Boresight  Angle  Rate 
3000  Feet,  30  Degrees. 
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Figure  19.  Best  Information  Horizontal  Trajectory, 
7000  Feet,  60  Degrees. 
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Figure  20.  Best  Information  Vertical  Trajectory, 
7000  Feet,  60  Degrees. 


Figure  21.  Best  Information  Boresight  Angle 
7000  Feet,  60  Degrees. 


Figure  22.  Best  Information  Boresight  Angle  Rate 
7000  Feet,  60  Degrees. 
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Figure  23.  Best  Information  A z  and  El  Angles 
7000  Feet,  60  Degrees. 
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Figure  24.  Best  Information  Az  and  El  Angle  Rates 


7000  Feet,  60  Degrees. 


Figure  30.  Minimum  Time  Az  and  El  Angle  Rates 
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